The evolution of obscured accretion 



Roberto Gilli*, A. Comastri*, C. Vignali f , P. Ranalli f and K. Iwasawa* 

*INAF - Osservatorio Astronomico di Bologna, Via Ranzani, 1 - 40127, Bologna, Italy 
1 Universita degli Studi di Bologna, Dip. di Astronomia, Via Ranzani, 1 - 40127, Bologna, Italy 

Abstract. 

Our current understanding of the evolution of obscured accretion onto supermassive black holes 
is reviewed. We consider the literature results on the relation between the fraction of moderately 
obscured, Compton-thin AGN and redshift, and discuss the biases which possibly affect the various 
measurements. Then, we discuss a number of methods - from ultradeep X-ray observations to the 
detection of high-ionization optical emission lines - to select the population of the most heavily 
obscured, Compton-thick AGN, whose cosmological evolution is basically unknown. The space 
density of heavily obscured AGN measured through different techniques is discussed and compared 
with the predictions by current synthesis models of the X-ray background. Preliminary results from 
the first half of the 3 Ms XMM observation of the Chandra Deep Field South (CDFS) are also 
presented. The prospects for population studies of heavily obscured AGN with future planned or 
proposed X-ray missions are finally discussed. 
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BH/GALAXY CO-EVOLUTION 

The observed scaling relations between the structural properties of massive galaxies - 
such as their bulge mass, luminosity, and stellar velocity dispersion - and the mass of 
the supermassive black holes (SMBHs) sitting at their centers suggest that the build-up 
of galaxies and the growth of black holes are closely related, which is often referred to 
as BH/galaxy co-evolution. Since the BH growth is thought to happen primarily through 
efficient accretion phases accompanied by the release of kinetic and radiative energy, 
part of which can be deposited into the galaxy interstellar medium, Active Galactic 
Nuclei (AGN) are believed to represent a key phase across a galaxy's lifetime. Support 
for this hypothesis comes from several lines of evidence such as i) the match between 
the mass function of SMBHs grown through AGN phases and the one observed in 
local galaxies [1, 2]; ii) the cosmological "downsizing" of both nuclear activity and 
star formation [3, 4]; iii) the detection of large quantities of gas outflowing from active 
nuclei (see [5] for a recent review). A number of semi-analytic models to explain the 
BH/galaxy co-evolution have been proposed over about the past decade [6]. These 
models follow the evolution and growth of dark matter structures across cosmic time, 
either through the Press-Schechter formalism or through N-body simulations, and use 
analytic recipes to treat the baryon physics within the dark matter halos. A common 
assumption of these models is that mergers between gas-rich galaxies trigger nuclear 
activity and star formation. Recently, a BH/galaxy evolutionary sequence associated to 
"wet" galaxy mergers has been proposed [7], in which an initial phase of vigorous star 
formation and obscured, possibly Eddington limited, accretion is followed by a phase in 



which the nucleus first gets rid of the obscuring gas shining as an unobscured QSO, then 
quenches star formation, and eventually fades, leaving a passively evolving galaxy. 

THE OBSCURED AGN FRACTION 

The various BH/galaxy co-evolution models have been successful in many respects, 
being able to explain the local BH/galaxy scaling relations, the local BH mass function, 
the QSO luminosity function and clustering [7, 8]. Further testing is however needed, 
given the large numbers of parameters involved in the modeling. A built-in feature of 
the BH/galaxy evolutionary sequence is that an obscured accretion phase preceeds a 
clean accretion phase, at least in powerful, QSO-like objects. Is it therefore reasonable 
to expect that the fraction of obscured AGN was higher in the past? This depends on 
many parameters such as i) the physical scale of the absorbing gas and how it is driven 
towards the BH; ii) the relative timescales of the obscured and unobscured phases; 
iii) whether the absorbed-to-unabsorbed AGN transition occurs also in low-mass/low- 
luminosity objects (i.e. Seyfert galaxies). 

Not many theoretical predictions on the evolution of the obscured AGN fraction 
are available in the literature. Depending on the assumptions, some models predict an 
increase of the obscured AGN fraction with redshift [9], some others do not [10]. 

From an observational point of view the situation is also debated. An increase of the 
obscured AGN fraction with redshift among X-ray selected AGN has been for instance 
observed by [11, 12, 13, 14], but other works did not find any evidence of this trend 
[3, 15, 16]. Much of the uncertainty can be related to the lack of statistically large 
samples of AGN with high-quality X-ray spectra. Indeed, using small samples which 
cannot be split into narrow luminosity and redshift bins, it is often complicated to 
disentangle any redshift-dependent from luminosity-dependent effects (it is generally 
agreed that the obscured AGN fraction decreases with luminosity, but the slope and 
normalization of this relation differ from paper to paper). Furthermore, the absorption 
column density is often estimated either from the X-ray hardness ratio, or from the 
absence of broad emission lines in the optical spectrum. Both methods suffer from 
caveats and limitations [17], and ideally the absorbing column should be measured 
through X-ray spectra with good photon statistics. 

An attempt to measure the evolution with redshift of the obscured fraction of AGN 
as a function of their X-ray luminosity has been recently done by Hasinger [13], who 
compiled a sample of ~ 1300 AGN selected in the 2-10 keV band and split it into 8 lumi- 
nosity bins. In Fig. 1 (left) is shown the evolution of the obscured AGN fraction at lumi- 
nosities around L2-10 ~ 10 44 erg s" 1 as measured in [13]. Although some approximated 
correction has been performed (see Section 6.2 in [13]), this is essentially an observed 
fraction which still suffers from selection biases against obscured objects: because of K- 
correction effects, this bias is less severe towards higher redshifts, when the photoelectric 
cut-off moves towards lower frequencies. We verified the strength of K-correction effects 
by considering the X-ray background (XRB) model by [16], which matches the obscured 
AGN fraction at logL2-io = 44, and in which this fraction is constant with redshift. By 
folding the logN-logS relations of obscured and unobscured AGN with logL2-io = 44 
with the sky coverage curve of [13] we found the trend shown in Fig. 1 (left): an intrin- 




Redshift redshift 



FIGURE 1. Left panel: fraction of obscured AGN with L2-10 ~ 10 44 erg s -1 as observed by [13] 
(Hasinger 2008, datapoints) compared with the model expectations by [16] (GCH07, solid lines). The 
intrinsic fraction of obscured AGN with L2-10 — 10 44 erg s -1 , assumed to be constant with redshift in [16], 
is shown as an horizontal line. When folding this intrinsic fraction with the sky coverage curve of the AGN 
sample in [13], because of selection and K-correction effects, the observed bent solid curve is obtained. 
Right panel, space density of CT AGN as measured through [O III]5007 selection (datapoints at z < 0.7) 
and through IR selection (z > 0.9). Model curves for CT AGN with different X-ray luminosity (as labeled) 
are from [16]: each datapoint should be compared with the curve marked by the same, smaller symbol. 
The shaded region shows the redshift range accessible to [Ne V]3426 selection in optical spectroscopic 
surveys; CT AGN candidates selected through X-ray spectroscopy in the XMM-CDFS cover the same 
redshift range and further extend to higher redshifts. Adapted from [26]. 

sic obscured fraction which is constant with redshift translates into an observed fraction 
which is increasing with redshift. Thus, the trend at logL2-io = 44 found in [13] might 
be mostly driven by K-correction effects rather than by a truly evolving obscured AGN 
fraction. The same conclusion also holds at lower luminosities. A possible way to recon- 
cile this result with the redshift evolution found in [1 1] and [12], who take into account 
selection effects in their estimates but do not split their samples into narrow luminosity 
bins, is to hypothesize that the obscured fraction increases with redshift only for very lu- 
minous QSOs. Indeed, the luminosity dependence of the obscured AGN fraction in [16] 
appears to be too shallow when compared with recent results [13, 17], overestimating the 
number of obscured QSOs with logL x ~ 45: a too large number of local obscured QSOs 
assumed in [16] might mask any increase with redshift. The possibility of a stronger 
evolution of the obscured AGN fraction for more luminous sources has, for instance, 
been proposed in the blast- wave semi-analytic model by [9]. 

THE MOST HEAVILY OBSCURED AGN 

If the overall evolution of obscured AGN is uncertain, the evolution of the most heavily 
obscured and elusive ones, the so-called Compton-Thick AGN (Nh > 10 24 cm -2 ; here- 
after CT AGN), is completely unknown. Only ~50 bona-fide CT AGN are known in the 



local Universe [18, 19], but their abundance has nonetheless been estimated to be com- 
parable to that of less obscured ones [20]. Also, synthesis models of the XRB predict 
that from ~10 to ~30% of the XRB peak emission at 30 keV is due to the emission of 
CT AGN integrated over all redshifts. In the absence of any information, the luminosity 
function and evolution of CT AGN in XRB synthesis models have been usually assumed 
to be equal to those of less obscured ones. 

In recent years there have been many attempts to constrain the space density of 
CT AGN in different luminosity and redshift intervals exploiting different selection 
techniques. Very hard (> 10 keV) X-ray surveys are still limited in sensitivity and are 
just sampling the local Universe [21]. As a result, the population of local CT AGN 
detected by INTEGRAL and Swift is producing only a tiny fraction (< 1%) of the 
XRB. To sample AGN at higher redshifts, and in particular at z ~ 1, where the bulk 
of the "missing" XRB is thought to be produced, one needs to rely on deep X-ray 
surveys in the 2-10 keV band and try to select CT AGN either directly from X-ray 
spectroscopy, or by comparing the measured, obscured X-ray emission (if any) with 
some other indicator of the intrinsic nuclear power: IR selection can track the nuclear 
emission as reprocessed by the dusty absorber; high-ionization, narrow optical emission 
lines can sample scales free from the nuclear obscuration. X-ray stacking of IR-selected, 
but X-ray undetected, sources has been used to estimate the space density of CT AGN 
at z ~ 1 - 2 [22, 23, 24, 25]. The comparison between the [O III] 5007 and X-ray flux 
has been used to select X-ray underluminous QSOs and then estimate the density of CT 
QSOs at z = 0.5 [26]. A summary of these and other measurements, compared with the 
predicted space densities of CT AGN by [16], is shown in Fig. 1 (right). To sample 
the population of CT around z ~ 1, in [27] we recently devised a selection method 
based on the [Ne V]3426 emission line, which can be applied to optical spectroscopic 
surveys with deep X-ray coverage to search for CT objects at z > 0.8, i.e. at redshifts not 
reachable with [O III]5007 selection. We calibrated an X-ray to [Ne V] flux diagnostic 
ratio (X/NeV) based on a sample of ~ 70 local Seyferts, verifying that objects with 
X/NeV<15 are almost invariably Compton-Thick. We then applied this diagnostic to 
several samples of obscured and unobscured QSOs: for the objects at z < 0.8 presented 
in [26] we found an excellent agreement between the absorption estimated from the 
X/OIII and the X/NeV ratio. Furthermore, despite the very low statistics, « 50% of the 
[Ne V]-selected type-2 QSOs at z > 0.8 appear to be good CT candidates. 

[Ne V], [O III], and IR selection are all indirect ways to select CT AGN, since the 
CT nature of an object is inferred from the faintness of its X-ray emission relative 
to an indicator of the intrinsic power, and then suffer from a number of systematics 
like reddening in optically selected sources and non-simultaneous data acquisition. To 
sample bona fide CT AGN over a broad redshift range and, in particular, at z ~ 1, very 
deep X-ray exposures are needed to collect photons and allow X-ray spectral analysis. 
This is indeed one of the goals of the deep XMM survey in the CDFS (P.I. A. Comastri). 
The exposure time, cleaned of background flares, currently amounts to ~ 1.8 Ms and 
by the end of the program, foreseen by March 2010, should reach ~ 2.5 Ms. A number 
of CT candidates, identified in the 1 Ms CDFS catalog on the bases of their flat (low 
quality) spectrum [28], are being confirmed as such by the higher quality XMM spectra 
(Comastri et al. in prep.), including the well-known CT candidate CDFS-202 at z=3.7 
[29] . Further significant improvements are expected by the forthcoming DDT extension 



of the Chandra CDFS exposure from 2 Ms to 4 Ms. 

FUTURE PROSPECTS 

Current deep (Ms) X-ray exposures are limited to very small areas (< 0.4 deg 2 ), and do 
not allow construction of sizeable samples of distant, bona-fide CT AGN. A number 
of X-ray missions are being developed or proposed for the near to mid-term future 
which should address this issue: NuSTAR, ASTRO-H (approved by NASA and JAXA, 
respectively), EXIST and NHXM (proposed to NASA and ASI, respectively; see these 
Proceedings) will primarily focus on the high-energy (> 10 keV) radiation, which is 
the least affected by absorption. Their sensitivity is expected to be up to ~ 2 orders of 
magnitude better than current instrumentation. The redshift distribution of CT objects 
collected by these missions will peak at z ~ 0.3 - 0.4, and will allow population studies 
of CT objects at redshifts z < 1. The proposed International X-ray Observatory (IXO, 
[30]), thanks to its large effective area (3 m 2 at 1 keV) has the potential to provide 
very high quality spectra for distant CT AGN, but the small (< 0.1 deg 2 ) field of view 
(FOV) might represent a limitation to collect large object samples. A mission which is 
certainly expected to allow population studies of distant (up to z ~ 4) CT AGN, thus 
providing a measure of their cosmological evolution, is the proposed Wide Field X-ray 
Telescope [31], which has a unique combination of wide FOV (1 deg 2 ), large effective 
area (1 m 2 at 1 keV), sharp angular resolution (~ 5" - 10") constant across the FOV, and 
low instrumental background to provide deep X-ray surveys in the 0.5-7 keV band over 
large portions of the sky. The equivalent of -2000 XBootes fields, ~ 3000 C-COSMOS 
fields and ~1000 2 Ms CDFS fields is expected to be returned by WFXT over 5 years 
of operation. As an example of these surveys' potential, based on the current design and 
observing strategy, the WFXT project is expected to return a sample of ~ 500 objects at 
z > 1 which are bona fide CT AGN, i.e. with more than 500 net counts in the 0.5-7 keV 
band 1 (the number of simple detections of CT AGN will be obviously much larger). 
In Fig. 2 we show two such objects (CDF-202 and CDFS- 153, another CT AGN, at 
z = 1.53, found in the XMM-CDFS) simulated using the WFXT response matrices for 
the goal design. It is evident that 500 X-ray photons are sufficient to unambiguously 
reveal their CT nature. Furthermore, WFXT is expected to return a global sample of 
300,000 AGN in the redshift range z = - 5 with even larger photon statistics (> 1000 
net counts), for which a robust measurement of the column density is possible. Among 
the many applications, such an object sample will definitely allow measurement of the 
evolution of obscured accretion as a function of AGN luminosity. 
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For comparision, depending on the observing strategy and background level, in 5 years of operation 
IXO is expected to return from <100 bona-fide CT AGN at z > 1 to a sample as large as the WFXT one. 



CDFS-202, z=3.700 „ CDFS-153, z=1.536 




FIGURE 2. 400ks WFXT simulated spectra of CDFS-202 and CDFS-153, two high-z CT AGN with 
/ 2 _io > 10" 15 erg cm -2 s _1 observed in the XMM-CDFS. The 100 deg 2 WFXT deep survey [31] is 
expected to reveal about 500 objects like these - in terms of obscuration and photon statistics - at z > 1- 
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